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Abstract. Rabi oscillations of two- level system driven by pulse train is a basic 
process involved in quantum computation. We give a full quantum treatment of this 
^ | process, and show that the population inversion of this process collapses exponentially, 

has no revival phenomenon, and has a dual-pulse structure in every period. As 
an application, we investigate the significance of above results in ion-trap quantum 
computation. We find that when the wavelength of the driving held is of the order 10 -6 
m, the lower bound of failure probability is no less than 10~ 2 after 10 2 operations. It 
qq | means that in Cirac-Zoller scheme, after about 10 2 CNOT gates, the failure probability 

0^ ' arrives at 10 -2 . This value is about the generally-accepted threshold in one error- 

. correction period in fault-tolerant quantum computation. 
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1. Introduction 

The presented quantum algorithms show that quantum computer (QC) can solve several 
famous problems intractable on classical computers [I], then challenge most public- 
key cryptosystems in use [2j [3] . Many proposals for implementing QC have been put 
forward. Among them, the cold ion trap scheme (Cirac-Zoller scheme) [I] is the earliest 
and most promising one, e.g., a scalable, multiplexed ion trap for quantum information 
processing has been demonstrated [5j. Implementation of quantum logic gates in this 
scheme is realized via Rabi oscillation of ions driven by pulse train of laser fields. The 
interaction of a single atom with radiation field is a basic interaction in physics. In [6], a 
nonperturbative fully quantum-theoretical analysis describing the transient spontaneous 
emission of an initially excited two-level atom in a one-dimensional cavity with output 
coupling is presented. In [7], they present observations of quantum dynamics of an 
isolated neutral atom stored in a magneto-optical trap. 

The theoretic measure adapted in pE] is a typical one which takes the laser field 
as a classical field. However, considering the quantum nature of the driving field, one 
may get results which is different from the result under classical treatment. There are 
generally two ways to take account of the quantum nature of the field. One is to add 
quantum fluctuations to the classical treatment j8] . However, there are many operations 
in quantum computation, then in this case whether this method is suitable is still to be 
considered. The other way is to quantize the field and work out the result [9] . To do this, 
we first consider the Rabi oscillation driven by quantized pulse train. This is a basic 
atom-photon interaction process, and quantum computation is only one application 
of it. Then we can analyze the failure probability in ion-trap QC, and have further 
discussions. 

Rabi oscillation driven by a quantized continuous-wave (cw) field, accompanied by 
collapse- revival phenomenon [TQJ [Til 1121 EE], is a typical phenomenon of atom-photon 
systems. However, Rabi oscillation driven by quantized pulse train has not been fully 
investigated. Conceivably, it may have different phenomena compared with that driven 
by a cw field. 

Fault-tolerant quantum computation (FTQC) allows the computer to work 
normally even when its elementary components are imperfect. However, the threshold 
theorem in FTQC requires the failure probability of each components below some 
threshold [H]. Then we can compare the failure probability with the value of threshold, 
and reach some meaningful conclusion in quantum computation [9]. 

This paper is arranged as follows: in Section [2} we describe a method to deal with 
quantum transformation of a two-level system after one coherent pulse, which expresses 
the relationship between the density matrices for the two-level system before and after 
one coherent pulse. In Section [3] we investigate the properties of Rabi oscillation driven 
by pulse train. In Section H] we study this kind of Rabi oscillation in ion-trap quantum 
computation, and get the failure probability. In Section |5] we give out some discussions. 
In Section [6] some conclusions are reached. 
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2. Quantum transformation of a two-level system involving one coherent 
pulse 

2.1. Modeling 

The two-level system driven by repeated pulses is an open system, the usual way to deal 
with such a system is Kraus summation and master equation method. However, for the 
specific problem here, which cannot be easily solved with those methods, we study it in 
this way: after a single pulse, we get the density matrix for the whole system (including 
a two- level system and the laser field), then obtain the reduced density matrix for the 
two-level system. We can get the relation for the state of the two-level system before 
and after the pulse, then the state of the two-level system after repeated pulses can be 
get. 

In [U [15], they use the Jaynes-Cummings model (JCM) [16] for the case where an 
atom in the free space interacts with laser field. However, the JCM is the model to 
describe the interaction of atom and single-mode field in a cavity. Actually, there are 
some discussions [TTt [T8| IT9] on the validation of the JCM in the multi-mode case. For 
example, in the paper by Enk and Kimble [TH], at section 2.3 "Atom-light interaction", 
they considered the case where an atom in free space interacts with a laser field, making 
use of Hamiltonian of JCM in Eq.(lO). They also pointed out that the Hamiltonian in 
Eq. (10) in the paper is valid for atoms in free space for less than one Rabi period, but 
there is no strict proof. 

We analyze the issue as follows: the sources of decoherence can generally lead 
to certain amount of failure probability on a single qubit or a pair of qubits. After 
many operations on the same qubit (or the same pair of qubits), the failure probability 
will generally accumulate to reach the threshold in the threshold theorem of FTQC. 
The corresponding operation number is the upper bound of operation number in one 
error-correction period when the given source of decoherence exists. Each source of 
decoherence considered can generally give corresponding upper bound of operation 
number, but this upper bound can be enlarged by technique improvement. When the 
upper bound increased to a huge value, it is actually not a bound. 

The decoherence caused by field quantization can also give its upper bound of 
operation number. It is obtained theoretically, thus cannot be enlarged by technique 
improvement. The calculation of this decoherence should include the interaction of all 
modes in the radiation field with the two-level system. When using JCM, only one mode 
of the field is considered, and can also give an upper bound of operation number 0x- 
The accurate upper bound of operation number from field quantization /3 < flx, because 
the spontaneous emission induced by vacuum modes is not considered in JCM. Then if 
we use JCM to estimate the upper bound of operation number in one error-correction 
period from field quantization, we can get meaningful results. The two-level system 
driven by pulse train can be described as 

H = hg (e^a+a + e _i V<r_) , (1) 
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where g is the coupling constant, is the beam phase, a + and <r_ are the raising and 
lowering operators of the two-level system, and a) and a the creation and annihilation 
operators of photons, respectively. Then the unitary time-evolution operation is given 
by 

U{t) = cos (gty/ala + l) |1>{1| + cos (gts/ata) |0)(0| 

sin [gtyfcfia + 1 ) sin ( gt\J a^a + 1 ) 

-a|l)(0| + e-'^t \ y |0)(l| 



V ata + 1 



(2) 



with |0) and |1) the ground and excited state of the two-level system respectively. 

Generally, the initial state of the whole system is \if) (0)) = J2^=o c n\ n ) ® 
(a|0) + where |c n | 2 = 6 " n , and |a| 2 + |/3| 2 = 1. A single qubit gate is usually 

implemented through a /c7r pulse in Cirac-Zoller scheme, whose duration t satisfies 
gt y/E = *y [T5], with n the mean number of photons in the pulse. After a kn pulse, 
the state for the two-level system and laser field is 




kir\/n.. , ;i .kitx/n^, 
( (*s(— l=r)\0,n) -ie^sin(— 1) 



2Jn 



+/?< J2 Cn 



. n=0 



COS 



|l,n) — le v sm( = — )|0,n+l) 



2y/E 



The corresponding density matrix for the state in (jHJ) is p^L; = 



• (3) 



This 



matrix contains the information of both the two-level system and the field, but we are 
interested only in the two-level system. Thus we obtain the reduced density matrix 

\a\ 2 S 4 + - a*(3)e^S 2 + |/3| 2 (1 - S 6 ) af3*S 5 + i(|a| 2 e i( ^i - \(3\ 2 e-'^S 7 ) + a*f3S 3 

a*pS 5 - i(H 2 e i( ^i - \/3\ 2 e-'^S 7 ) + a/3*S 3 \a\ 2 (l - S 4 ) - ±(a/3* - a*/3)e i *S 2 + \/3\ 2 S 6 



P 



(i) 



n=0 



n! V w + 1 



n=0 

oo 

s 3 = E 



e n 



ft! 



n=0 



s 4 = E 

n=0 

oo 

S = = E 



e n 



n kix^fn kTTy/n + 1 , 
cos( — ) smi 



kn 



2(n + l) 



sin 



2v^ 2^ 
,kn\/n + 1. 



n n n / n .k^Jn .kny/n + 1. 

— I — \/ 7 sm ( 7^ J sm ( /= — 

n\ V n+ 1 2^ 2^ ' 



cos 



n! 



,kny/n 

~2~Vn 



n=0 



s. = E 



n=0 



nl 



rv. 



k-nJn kny/n + 1, 

cos! =r- ) cos( = — 

v 2^n ' y 2^n ' 

. kiry/n + 1. 



cos 



2^fi 



-)• 



(4) 
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n=0 



n! 



n .kizyn + 1. .kix^/n 

- cos ( — 7^ ) sm (^-7^)- 

n 2\m 2vn 



Transforms of the density matrix after a coherent pulse 

Consider the relationship between p^ and the density matrix of corresponding initial 
state p(°) = \ip (0))(ip (0) \. For a two-level system, the density matrix p satisfies 
the condition p = ~(I + r • c) [2], r is the Bloch vector for state p, \r\< 1, 



a 



a x a y a z 



Let r( r 



(m) (m) (m) 



denotes the Bloch vector of p^ m \ An arbitrary 



trace-preserving quantum operation is equivalent to a map of the form r — > r' = 
Mr + c [H], here M and c contain the properties of the system and are independent 
of the state. Based on this, it can be seen that = Mr^ + c, here c = 



Sre-* - SW* S 4 + S 6 -l 
M = 

then r( m ) = Mr {rn - 1] + c 



s 3 + s 5 







Mi 



Mi 



5 5 -5 3 -(e^ + e-^r) 
5 2 e^ S 4 - 5 6 



^.5. Calculation of the sums in the density matrix 

It is necessary to get accurate values of Si (i = 1, • • • , 7) to evaluate the behavior of pulse 
train. The usual algorithm (saddle-point approximation) can only reach a precision of 
1 / \ffi. Our algorithm achieving any given precision instead of the usual algorithm is as 
follows. 

Suppose n is not small, for the sum 

Y] , fio(n,n,k), 
^— ' n! 

n=0 

(1) Substitute n in fio(n, n, k) with (x+l)n , we get fn(x, n, k) = f i0 ^(x+l)n, n, kj . 

(2) Do Taylor expansion to x p for fn(x, n, k) at x = 0, and get fefa, n, k). 

(3) Since sum Yl^Lo 6 "' ) " , " n fc can be get accurately, we replace x in fi2(x,n, k) by 
^ and get f a (n,n,k). 

(4) Use f i3 (n,n,k) instead of fio(n,n,k) in the expression of Si(n,k) to calculate 
the new sum and get /^(n, fc). 

(5) Substituting n into / i4 (n, fc), we obtain a high-precision result of the original 
sum Si(n, k). The value for Si (i — 1, • • • , 7) in the cases where we expand fa(x, n, k) 
to a; 10 and x 15 are compared in Table HJ 

The precision of the sums Si(i = 1, 2, • • • , 7) is ensured by the following theorem: 
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Table 1: Values for Si(i = 1, 2, • • • , 7) for n = 10 4 and k = 2. "Valuer" denotes value 
of the resulting sums of the algorithm when we expand fn (x, n, k) to x 10 and "Value2" 
denotes that when we expand fn(x, n, k) to x 15 . Valuel and Value2 are the same to the 
precision 10 -23 . 



Sum 


Valuel 


Valuc2 


Si 


0.000 039 303 916 656 063 668 561 519 091 


0.000 039 303 916 656 063 668 561 194 770 


s 2 


0.000 039 265 164 255 300 772 996 074 590 


0.000 039 265 164 255 300 772 995 750 283 


s 3 


0.000 246 659 192 761 352 167 541 307 293 


0.000 246 659 192 761 352 167 542 402 758 


s 4 


0.999 753 309 972 685 637 856 777 333 369 


0.999 753 309 972 685 637 856 776 237 858 


s 5 


0.999 753 316 133 881 571 308 212 070 145 


0.999 753 316 133 881 571 308 210 974 684 


s 6 


0.999 753 322 301 165 250 291 025 614 276 


0.999 753 322 301 165 250 291 024 518 866 


s 7 


0.000 039 226 416 698 193 975 826 600 887 


0.000 039 226 416 698 193 975 830 095 264 



Theorem 1: For every given integer I << n, let 
In (y2n l -^(l + l)lnn) 



P 



s In n — In ( (I + 1) Inn 



(5) 



1 (Z + l)lnn /(/ + l) 2 (lnn) 2 

«o = -7= + LJL ^= + \ + 2(1 + 1) Inn. 

\/n \/n V n 



If a < a « \/n, then 

yv e — _^ _ e — n _^ + / 1 

n=0 n=0 

here p, fio(n,n, k), fi^(n,n, k) are parameters defined in the algorithm above. 



(6) 



Theorem 1 can be proved using the following three lemmas (see Appendix A for 
the detailed proof): 

Lemma 1: For every given a « y/E, 



X] „, (fio(n,n,k) - f i3 (n,n,k)) = o 



711 



p+1 



(Vn) 



p-i 



(7) 



n=n— ay/n 

here p, f i0 (n,n, k), f i3 (n,n, k) are parameters defined in the algorithm above. 



Lemma 2: For every given integer I « n, a « \/n, if a > \J (I + 1) Inn, then 

1 



E 

n=0 



e -r < —i 

n\ rv 



(8) 



where k = |"n + o;\/n] . 

Lemma 3: For every given integer I << n, a << \/E, if 



1 (Z + l)lnn /(/ + l) 2 (lnn) 2 nn iN1 

n > ^= + - t4 + \ ~ — - + 2(l + l)\nn, 

'n \Jn V n 



(9) 



then 
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where k' — [n + oja/^J . 

We expand f\(x,n, k) to a; 15 at x — 0, and find the value of the sum is the same 
to the precision 10~ 23 as that when we expand fi(x,n,k) to x 10 . It implies that the 
precision of the sum is much higher than Eq. ([5]) shows. The reason is probably that we 
have not considered of the periodicity of trigonometric functions, and the precision of 
the sum may be remarkably improved since the positive and negative terms will cancel 
each other out. 

For small n, we need only to require t satisfying 

it- 1)! > e- fl n t+l , (11) 

where t is the parameter in sum Si(n, k) = Yln=o £_ ^f~^°( n ' ^> ^) + °\Wj ■ ^ or a §i ven 
precision I, we can search out the smallest t satisfying (11 II) . e.g., when n = 10 and 
I = 20, we set t = 55. 



3. Population inversion 



3.1. Final state of the two-level system after pulse train 
Provided r^ = Mr {m ~ 1] + c , then 

r (m) = M m r (0) + ( M m-l _| + M + I)c. 

It can be seen from Sec. 12.21 that 



M r - 



For any real matrix A 



(S 3 + S 5 ) m 
O M™ 



S 5 -S 3 -(St + Sr) 

S2 5*4 + Sq — 1 



a b 
c d 



we obtain (see Appendix B) 



Am ^ r A - 

A m = —^—1 + — 



(m) 



2ig 



-X 26 
2c K 



where A 



(m) 



(12) 



(13) 



± 



A™ ± A™ with Ai and A2 eigenvalues of A, K = d — a,Q 



-i\/(a - d) 2 + 46c. When (a - d) 2 + Abe < (which is the case for Mi) 

J 



A m = \X\ m cos(m6)I + sin(m6 



VdetJ 



where |A| 2 = ad - be, sin9 = \ J2 - a2 +£+f c , J 



m—l 



a — d 2b 
2c d — a 



. Therefore, 



m—l 



I + Mi + ■ ■ ■ + M™- 1 = \J2\X\ j cos(j9) I + [ \ X \ j sinO'0) 



J 



j=0 



j=0 



VdetJ 
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Since 



we have 



thus 



then 



m— 1 



m—1 



m—1 



Y,W\«*W +ism(j0)]} = £(|A|V^) = £(1*1**)' 

1 - |A|e- ie - \\\ m e< m °) + \\\m+i e i( m -i)e 
~ 1 + |A| 2 — 2|A| COS6* ' 

1 - |A|cos#- |A| m cos(m#) + |A| m+1 cos(m - 1)9 
1 + |A| 2 -2|A|cos# ' 

|A|sin0- |A| m sin(m#) + | A| m+1 sin(m - 1)9 
1 + |A| 2 -2|A|cos# ' 



m—1 

£|A|'"cos(j0) 

3=0 
m—1 

3=0 



I + M 1 + --- + M r 1 
1 



m—1 



(1 - |A|cosfl 



1 + |A| 2 — 2|A| cos6» 
+ (|A|sin0- |A| m sin(m0) + |A| m+1 sin(m - 1)6) 

= B[ m) I + Bi m) J, 



cos(m6) + |A| m+1 cos(m - 1)6)1 
J 



VdetJ 



r O) — 



+ 



(s 3 + s 5 y 
o 

l-(S 3 +S 5 ) m 
1-S 3 -S B 



cos(m6)I + sm(m6)^=L= 



,(o) 



O B^I + B^J 



c. 



3.2. Population inversion after pulse train 

Suppose the initial state is |1), if we have applied kn pulses for m times, the population 
inversion is 

W m = \(1 - r^) - l -a + r<T>) = -r<r\ 



we have 
W m = |A 



cos(m6) + sin(m#)- 



322 



v^etjJ 1+ |A| 2 -2|A|cos£ 



x{ [|A| sind - |A| m sin(m0) + |Ar +1 sin(m - l)d\ / 21 (S 7 - S ± ) 
L -J v det J 

+ [(l - |A|cos0- |A| m cos(m0) + |A| m+1 cos(m - 1)^ 



+ (|A| sin9 - |A| m sin(m0) + |A| m+1 sin(m - l)fl) -jt= = ] (S 4 - 5 6 )}. (14) 
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To get the inversion between the mth and (m + l)th /c7r pulse, we should first get 
the corresponding density matrix of the two-level system 

Pm (t) = tr {f/(t)[p (m) ® Pi]U\t)} , (15) 

where U(t) is the unitary time-evolution operator mentioned earlier, and p\ is the density 
matrix for the laser field. A detailed calculation gives the probability that the ion is in 
state |0): 

where S, ^ = £~=o e -^f cos 2 gty/E, S 9 = ^ gtyfr+T and S l0 = 

Yl^=o 6 rX y/n s ^ n ^Qfcy/n, which can all be got with high precision using our algorithm 
in Section [231 and the inversion is W m (t) = 1 — 2p(t), < t < 

Now think of the difference of the oscillation driven by pulse train and a cw field. 
The population inversion for repeated 2tt pulses is shown in Figure [la] (given n = 10). 
We find there is a dual-pulse structure in every period, where the amplitude starts to 
increase from the point of 2tt pulses. The inversion decreases exponentially, different 
from a Gaussian function collapse envelope driven by a cw field. Besides, there is no 
revival phenomenon, but a small nonzero amplitude exists. The reason for this behavior 
is probably that the system we consider here is an open quantum system with dissipation. 

The inversion at the points of 2ir pulses when n = 10 4 is plotted in Figure El 
Results of fitting is 1.0031e- a0002ArR , 1.0193e- a0003ArR , 1.025e- a0005ArR for k = 1/2,1,2 
respectively, here is the number of Rabi periods. 



Pit) = « (Ss + S 9 ) + ri m \t)(S 8 - S 9 ) + rW(t)S w 



4. Failure probability of gate operation realized through Rabi oscillation 
driven by repeated pulses 

4-1. Estimation of n 

n which determines Si is an important parameter in our discussion. To estimate the 

mean number of photons in one pulse, we assume a fictious pulse propagating in 

the opposite direction at the same time. They form a standing wave instantly while 

overlapping in space. It can be seen that the mean number of photons in each pulse is 

about half of that in the standing wave. We now focus on mean number of photons in 

the imagined standing wave. 

The electric field E can be expressed as E = £y/h~. £ is usually given as £=\J ^7 

[20] . where u is the frequency of the single mode in a cavity, and V is the volume of the 

cavity. It can be seen that V ~ Act, with A the cross-sectional area of the beam, thus 

n = ^j^E 2 . For a kir pulse, gt\ffi = g ~ j- = -0=, with p ~ eao the electric dipole 

moment of the ion, e the charge of a electron, and Bohr radius, then we get t = 

Then we have n = ^^E. 
, p 

When the laser beam interacts with an ion, not all the photon in the beam are 
effective. Think of the effective mean number of photons. When a laser beam is 



10 



1 



0.5 




V 
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x 

(b) 

Figure 1: Population inversion driven by different fields, given n — 10, r = gt. (a) 27r 
pulse train case. There is a dual-pulse structure in every period, where the amplitude 
starts to increase from the point of 2ir pulses. The inversion decreases exponentially, 
different from a Gaussian function collapse envelope driven by a cw field. Besides, there 
is no revival phenomenon, but a small nonzero amplitude exists, (b) Corresponding 
inversion driven by a cw field. 

applied to a trapped ion, the total resonant scattering cross section for an atomic dipole 
transition is a = 3A 2 /27r |21j . and the cross section for scattering out of the paraxial 
modes is a e ff = 3\ 2 /8n [22]. Then the effective interaction area is & e ff, and the photons 
in volume cr e ffct is effective. For each photon, the probability to be in area cr e // is 
and the probability are independent for the photons. It can be seen that the effective 
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k=1/2,l=4 

--- k=1,l=4 
k=2,l=4 




0.2 



1000 2000 3000 4000 5000 6000 7000 

Figure 2: Inversion at the points of 2tc pulses E(W) versus number of Rabi periods 
N R , n = 10 4 . Fitting results are 1.0031e- a0002AfR , 1.0193e-°- 0003iVR , 1.025e- a0005JVR for 
k = 1/2, 1, 2 respectively. 



mean number of photons is 

_Cr e ff ktO&eff^rn 

n e ff = n — J -± = 1J — E. (16) 

A A p ( J 

One case people are interested in is the sideband transition, where the laser detuning 

A = ±u; t , here u t is the frequency of the trap. Because of AC-Stark shift and off- resonant 

transitions, the sideband Rabi frequency Q + has upper bound [23]. People have adopted 

methods to partially cancel the effect, and it seems feasible to have Q + < u t for special 



temporal and spectral arrangements of the laser field [24] . Since fi + = jy 2MujS *"> w ^ n 
M the mass for a single ion, we have 



A 2M 3 



From [25] and [26], it can be seen that 
cclqE pE 



Q - 



e 2 



Ane Mz^ ' 

where z s is the order of the separation between ions and is typically 10 to 100 /iin. 
Suppose z s = £A, from ( TITl) and ( fT8l) . we get 



2y/2h, e 



2 



5 



E<— — ( )?M^f?A"5. (19) 

pn Aneo 



Substitute back to ( fTBI) . we get 



32a2vrx V e 

= 6 x 10 7 fcM-^-fA3. (20) 
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In the cases we consider, it is suitable to limit k < 2, 9u < M < 200u (u= 
1.66057 x 1(T 27 kg). For M = 9u, k = 2, we get 

n = 3.4 x 10 14 r^i 

We can see that a large A and a small £ result in a large n. The curves of lg(n) is plotted 
in Figure |3] versus parameter £ from 2 to 100. When A = 10~ 6 m and £ = 2, we get 
n = 2.3 x 10 3 . 




Figure 3: Logarithm of mean number of photons lg(n) as a function of £ and A. It can 
be seen that n increases with A and decreases with £. 

There are also authors calculating n in a kir pulse in another way [15] : they 
considered the situation where a single laser is used to drive Rabi oscillation of the 
atom, and adopted the formalism introduced by Blow et al. [27], taking the laser 
as a continuous-mode coherent state. They worked out the interaction time t for kir 
pulse as t = where d is the coupling constant of the atom and laser, and 

P is the power of the laser. Thus, the mean number of photons in one kn pulse is 

™ ~ /£z = utji\J wnere is the frequency of the representative single-mode 
coherent state. Thus, obviously, they take all the photons in area A as effective photons 
when considering the interaction, but actually the number of effective photons is much 
smaller. 

4-2. Accuracy of gate operation 

Suppose we have applied m times of coherent pulses and reached a state = 
i(7 + r^" 1 ) ■ a). Let |^) = a|0) + /3|1) be the expected state, the accuracy rate of 
gate operation realized through Rabi oscillation is 

p^ = ( a *(0\+P*(l\)p( m \a\0) + (3\l)) 

= i (\ _l r (o) r M + r (o) r M , (0) {m)\ 
o \ y y J 

= i(l + r (o). r M) ; (21) 



for a mixed state, \r^\ < 1, then p s < 1. The failure probability is 



Pi 



(m) 



P 
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detailed calculation results (see Appendix C ) 



Pi 



(m) 



-(r 



(0) . „(m) 



(ri 0) ) 2 [(5 3 + S 5 ) m - 1] + ((rf) 2 + (rl 0) ) 2 )[|A|^ cosM) - 1] 

+ |A|^ sin(m#)(det J)"* [((rf) 2 - (rf) 2 )(a - d) + rfVf (6 + c)] 

+ i?i m) (rj 0) c y + r(°) Cz ) + 4 m) [(r( )c y - rfc z )(a - d) + r<«V + r y °) Cz 6] 

Then we average over all initial states of the ion, and get the average failure 
probability. The failure probability for kn pulses with different n are shown in Figure HJ 
It can be seen that the failure probability increases with the number of Rabi periods 
A^ R and the value k, and is inversely proportional to n. 



(a) 



(b) 



k=2,l=4 

- k=1,l=4 
k=1/2,l=4 




100 




80 100 



Figure 4: Failure probability pf versus number of Rabi periods A^r, kn pulses are applied, 
and n is 10 4 (10 6 ) in a (b). The failure probability pf increases with the number of Rabi 
periods Nr and the value k, and is inversely proportional to n. 



5. Discussions 

5.1. The reason for tracing out the laser field after each pulse 

For the Cirac-Zoller CNOT scheme, there are five steps to implement it. The laser 
phase in the five steps are 7r/2, 0, 0, and — vr/2, which are correlated. However, there 
is no correlation between two different CNOT gates. Then after completing one CNOT 
gate, it is valid to trace out the laser field. The calculation in Sec. [2] for different pulses 
is an analog of different CNOT gates. This treatment is meaningful for our purpose, 
i.e., to estimate the lower bound of failure probability, or the upper bound of operation 
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number, after many CNOT operations. More-detailed analysis of the lower bound of 
failure probability after Cirac-Zoller CNOT operations can be seen in [28], and the main 



idea of [28J is summarized in Appendix D 



5.2. The permitted depth of quantum logical operation 

Failure probability we have calculated for the 7r sideband transition amounts to 10 -2 
after approximately 10 2 operations when n = 10 4 . For controlled-NOT (CNOT) gates, 
there are five steps in Cirac-Zoller scheme, and two steps are realized via Rabi oscillations 
driven by n pulses. Generally speaking, the failure probability after five steps is not 
less than that after one ir pulse. Then the failure probability after repeated 7r pulses 
is a lower bound of the failure probability after repeated Cirac-Zoller's CNOT gate. 
Thus the lower bound of the failure probability is 10 -2 after approximately 10 2 CNOT 
operations when n = 10 4 . 

The threshold theorem in quantum computation declares that an arbitrarily long 
quantum computation can be performed reliably if the failure probability of each 
quantum gate is less than a critical value. Knill used numerical calculations and obtained 
a threshold of the order 10~ 2 [29] based on a fault-tolerant structure suggested by 
himself. P. Aliferis et al. reached a threshold of 1CT 3 with provable constructions [30] . 

A parameter called permitted depth of logical operation describing the property 
of a physical realization scheme of QC has been given [9]: considering that different 
number state components of the driving field lead to different oscillation amplitudes, 
which become uncorrelated gradually, we can see that failure probability of quantum 
logic gates has a theoretical limitation. Combing this limitation given by the quantum 
nature of the field with the threshold theorem in FTQC, we can get permitted depth 
of logical operation. This parameter limits the number of operations on any physical 
qubit in one error-correction period. Then the permitted depth of logical operation here 
is less than 10 2 . 



5. 3. Others ' proposals which may have different results 

For Rabi oscillation driven by microwaves, the failure probability may be much smaller 
because of a large mean number of photons, but it becomes difficult to individually 
address each ions. Although an additional magnetic field gradient applied to an 
electrodynamic trap may individually shift ionic qubit resonances [31], thus making 
them distinguishable in frequency space, whether it can improve the permitted depth 
of logical operation needs further investigation. 

There exists a two-qubit gate scheme totally different from the Cirac-Zoller gate, 
such as the scheme implemented by the NIST group [32]. In the scheme off- resonant 
excitations of the stronger carrier transition is absent, and this allows a greater gate 
speed thus a higher laser intensity. Besides, additional Stark shifts can be efficiently 
suppressed by choosing almost perpendicular and linear polarizations for the laser beams 
[33] . Hence, studies on this type of gate may lead to different results. 
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6. Conclusions 



Firstly, we investigate Rabi oscillation of a two-level system driven by pulse train. We 
develop an algorithm to solve the infinite summation, with a high precision which has 
never been reached by former authors. We find that in this kind of Rabi oscillation, there 
is a dual-pulse structure in every period. The envelope of population inversion collapses 
exponentially, different from a Gaussian function collapse envelope driven by a cw field. 
Besides, there is no revival phenomenon, but an approaching to a tiny platform. 

Secondly, we consider the application to gate operation in ion trap quantum 
computation. We give a lower bound of failure probability. Our result is: that when 
the wavelength of the driving field is of the order 10~ 6 m, the mean number of photons 
cannot be greater than 10 4 . Then, after about 10 2 CNOT gates in Cirac-Zoller scheme, 
the lower bound of failure probability is of the order 10 -2 . 
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Appendix A. PROOF OF PRECISION OF THE ALGORITHM IN SEC. 

Proof of Lemma 1: For n — a\M < n < n + aJfi, i.e. — % < x < -7=, after 
expanding fn(x, n, k) at x — 0, we get the result fi%(x, n, k) satisfying f i2 (x,n,k) = 
fio(n,n,k) + o(x p ). It can be seen that fi2(x,n,k) = f i3 (n,n,k), thus we have 
fio(n, n, k) = fz{n, n, k) + o(x p ), then Eq. (j7j) is proved □ 

Proof of Lemma 2: For every given n satisfying n < k + 1 < n, we have 

< (fc+l)!> tilUS 



k _„• k 



f.k+1 ^k+l 



n=0 n=0 v ' 

k\k. 



From Stirling's formula k\ = y2itk(~) ei2fc,0 < 6 < 1, we have 

e ~"^TT < (l) k e- n n k+1 = e fc ""n(-)- fc . (A.2) 
Substitute k in formula ( 1A.2I) with n — a\/E, we have 



n Jn 



When a > y^T+T) Inn, we have ne ° 2 < ^, inequality ()8]) is proved. □ 
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Proof of Lemma 3: It can be seen that 

= n =ik' 00 ^n-k' 



E 



-71 



n" 



~.— n 



n'* n" k'\ 



n! k'\ in 

n=k' n=k' 



E 



< e 



,n K 



k'\ 



■T(— 
^ K k' + i 

n=k' 



,71-fc' 



n k hf + 1 
~Mk' + l-n 



when k' > n + 4, i.e, , , fc t 1 _ < k' we have 

n' ' k'+l—n ' 
= fc'+l 



< { z ) k'e- n n k ' +l < ne~ n 
k'\ k 



k'-i 



(k'-iy: 



with Stirling's formula we get 



zk'-l 



e "k' -n-k'+l ^ ^„-n n ^ ^„-ii-k'-l _ ^„-n-k'-l 



(_Y e~ n fi k +i < fie~ n - 

V 6 n <n (k'-iy. 



< ne n 



ne n 



Let A = < 1, rj — n^ , we then have 



1 



ne 



ne n n k ' 1 = ne n n k ' 1 < -r ( 

i( k' — 1 



k'-i 



< 



77/ 



<£> (-) A -eA > 



n i+i 



<^> A(l - In 77) > 1 + In A. 



Let A = 1-A, withO < A < 1, from ln(l+x) < x-\x 2 (x < 0), we get In A < -A-|A 2 , 
then a sufficient condition of (^) A — eA > is: 

(l-A)(l-lnr ? )>l-A-iA 2 , 

which results in A > A , here A = — In 77 + \J (hiT?) 2 + 2 In 77. Let = n + a y^n~, 
we get 



\/n 
1 



+ 



1- A 

l + l 

n 



- 1 n+ 1 



Inn + 



— = — Inn + 2 — = — Inn | yn 

n I n 



by using = 1 + m 77 + ^/ (In?]) 2 + 2 In 77 and 77 = (n) » . Because A > A ct > a , 
we get a sufficient condition of Lemma 3: 



a > —7= + 



Z + l 



Inn + 




1 (Z + l) Inn /(/ + l) 2 (lnn) 2 n/1 

= -= + - £ + \/ - — + 2(1 + 1) Inn. 

Vn vn V n 

then Lemma 3 follows. □ 

Proof of Theorem 1: From Lemma 1,2 and 3 we get: for every given I << n, 
if a satisfies 
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1 (7 + 1) Inn 

+ ^ — j= — + 



(A.3) 



'(/ + l) 2 (\nn) 2 / _ 
+ \/- — -+2(Z + l)lnft < a « Vn, 



(J7D, © and (DUD hold. Then 



nl 



ifio{n,n) - /a(n,n)) 



n=0 

n—ayK _ 
n=0 

— ^2 — fl 

^2 n \ (/»o(n, ft) - /i3(w, ft)) + ^2 

n=fi J rOt\ffi 



nl 



n=0 



nl 



-fa{n) 



+ £ 



ifio(n,n) - fa(n,n, k)) 



n=n—a\j n 
oo 



Ee""n" , , , , / 1 \ / a p+1 
— f a n, n, k + o — ) + o ^ 

n=0 v ' 

Leto(i r )=o(^^ T ),weget 

In [\/2n l ~^(l + l)lnft 



P 



then S, 

ioo e -"n n 



~ In n — In y(l + 1) In n 

o(^t). Since we can get exact result of 
i 



12n=0 n\ fi3( n i n 



T,n=o 5_ sr-/«(n,w), we get Si with precision of ^J. □ 

Appendix B. CALCULATION OF Mi m 

Let 



where 



then 



M 



Mi 



M r ' 



S 3 + S 5 
O Mi 



(B.l) 



£5 — 5*3 — (Si + St) 


A 


a 6 


S2 5*4 + Sq — 1 




c d 



(S 3 + S 5 ) m 
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Figure Bl: A(r) = (a — d) 2 +Abc < versus r = gt, where t is the pulse width. Different 
A(t) results in different behavior of Rabi oscillation driven by pulse train. For the cases 
we consider, r < 1, we can see A(r) < 0. 



Let (1, X2i) T and (1, X22) T be the eigenvectors of Mi with corresponding eigenvalues 
Ai and A2, then 

X21 = ^r[d -a + V '(a - d) 2 + 46c], 



2b 1 



£22 = 777 [d - a - yf{a - d) 2 + Abe] , 



2b 

Xi = ^[a + d+ (a - d) 2 + Abe], 

A 2 = -[a + d- ^ {a - d) 2 + Abe], 

we have plot A(r) = (a — d) 2 + 46c versus r = gt in Figure 
the cases we are interested in (r < 1). 
Let 



A(t) is below zero for 



1 1 

x 2 i ^22 



thus 



then, 



Mi = T\ 



Ai 
A 2 



1 l 5 



Mf = Ti 



AT 1 
A^ 



-1 



X21 — ^22 



— ^22^l + ^2lA 2 Aj — A 2 
^21^22(-A™ + A") ^21 A^ - Z^A? 



Denote A™ ± A™ = Ai m) , d-a = K, ^(a - cf) 2 + 46c = iQ, we have 

Mr 



A M A ( m ) 

/-..'« _ A + _J _l_ A - 



-tf 26 
2c K 



2 " ' 2iQ 

It can be seen that |Ai| = | A2 1 , let Aj = \\\e ld , A2 = |A|e~ ie , using |A| 2 = A1A2 = ad — be, 
we get 



A^ m) = 2(ad-6c) m/2 cos(m#), 
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A M = 2i(ad - bc) m/2 sin(m#) 
where 6 satisfies sin# = 



2ad-Abc-a 2 -d? 



4(ad—bc) 



then 



M™ = |A| 



cos(m6)I + sin(m6 



J 



VdetJ 



where 



J 



a-d 2b 
2c d — a 



(B.2) 



Appendix C. CALCULATION OF r( m ) 

It can be seen from flT2), (JED, flRll) . flR2|) that 

r (0) . r (m) = ( r (0)jT M m r (0) + ( r (0)\T 



Using 



and 



we get 



(s 3 + s 5 r o 

O M[ 



EZo (s 3 + s 5 r o 
o Y,Zo m} 

l-(S 3 +S 5 ) 



r (0) + f r (0)\T 



l-(5 3 +S 5 ) 



O B[ m) I + Bi m) J 



(r^) 2 (S 3 + S B ) m + 5| m) (rf c y + rf c z ) + ((rf) 2 + (r^) 2 )|A|^ cos(m0) 



+ IAI 



smlmfc 



VdetJ 
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,(o) ro) 
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.(0) 
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.(0) 



• (CI) 



( r (o) Cy _ r (o) Cz )( a _ d) + rf c y c + rf c z 6, 



,(o) 
y 

,(o) 



((r (0) )2 _ (r (0) )2)(a _ rf)+r (0) r (0) (6 + c); 



r 



(0) . _ /„ (0)\2 



rf ) 2 [(S 3 + 5 5 ) m - 1] + ((rf) 2 + (r(°)) 2 )[|A|f cos(mtf) - 1] 
+|A|? sin(m#)(det J) "3 [((r y 0) ) 2 - (rf ) 2 )(a - rf) + rfrf (fc + c)] 
+B[ m \r^c y + rf c z ) + ^[(rfc, - rf c z )(a - rf) + rJV + r^c z b] + 1. 



Appendix D. More-detailed calculation for CNOT gate [28J 

Now, we consider the full-quantized treatment of Cirac-Zoller CNOT gate. The main 
idea is as follows: in the Cirac-Zoller scheme, the qubits are trapped ions and the phonon 
of the ion string's collective motion. For the CNOT gate, suppose xth ion is the control 
qubit, and the yth ion is the target qubit. Then the initial state l^ )) of the qubits 
(including xth ion, yth ion and the phonon) is a superposition of 12 computational basis 
states. Let the coefficients for the computational basis states be onii = 1, ■ • • , 12). For 
each step of CNOT gate, full-quantized treatment is applied, and the laser field is not 
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traced out. Finally we get the state after one CNOT gate. Then we trace out the laser 
field, and get the density matrix of the ions and phonon after one CNOT gate. It has 
144 elements, and the simplest element is 

p\ x = a. x [aib 2 o 4: — a 3 b 3 b 4 + a 3 b i b 1 + aibi&i — a 3 b 6 b w + a 3 b 7 b4 + aii> 2 04 + ai5 5 5ioJ 

+«3 [a 3 S w Si + a%SiSi — ctsSgS^ — axSsS^ — aiS 6 S w + S%(—a 3 Sg + ai<SV) + c^S'sS'io] 

where Si(i = 1, • ■ • , 10) are sums with the form 5^^1 e ,"" /(n, n), with n the mean 
number of photons in the pulse. The sums can be get using the algorithm in 12.31 

Consider the state transformation for the qubits after one CNOT operation with 
initial state p^ = |^^)(^^|. For the first step of CNOT gate, suppose the state for 
the qubits after this step is pa , the operator-sum representation of hyperoperator in 
quantum open system has the results p^ = S(p^) = ^2 i E i p^Ej. If we let 

M 1 = Ei ® (E\) T + E 2 ® (E\) t + ••- = #! ® E{+ E 2 ® E* + -- ■ , 

then p^ = Mip (0) . Here, p is straightening of matrix p. Similarly, if we let the 
corresponding matrix of the second to fifth step be \ ;..\ /.-,. then 

piV = Mp(°\ (D.l) 

here M = M 5 M 4 M 3 M 2 Mi . 

According to the previous calculation, p(°\ p^ are 12 x 12 matrices. Then 
p^ 1 ) are 144 x 1 column vectors. From Eq. (ID . 1 1) . we can see that the matrix M 
corresponding to one CNOT operation is of the size 144 x 144. We need to determine 
the 144 x 144 elements. From Eq. ( ID. II) . we can only get 144 equations. We have 
found other relation to get 144 x 144 equations, and finally get the matrix. Once M 
is obtained, from = Mp(°\ we can get the state after t times of CNOT operation 
p(*) = MpQ-V = ■■■ = M'p (0) , then p® is get. 

is a state with parameter oti(i = 1, • • • , 12), then for different initial state, the 
actual is different. We can get the failure probability of CNOT operations with the 
expected state and the actual state. The final results we reach is: the failure probability 
is of the order of 10~ 2 after 10 2 CNOT operations. 
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1. Introduction 

The presented quantum algorithms show that quantum computer (QC) can solve several 
famous problems intractable on classical computers [I], then challenge most public- 
key cryptosystems in use [2j [3] . Many proposals for implementing QC have been put 
forward. Among them, the cold ion trap scheme (Cirac-Zoller scheme) [I] is the earliest 
and most promising one, e.g., a scalable, multiplexed ion trap for quantum information 
processing has been demonstrated [5j. Implementation of quantum logic gates in this 
scheme is realized via Rabi oscillation of ions driven by pulse train of laser fields. The 
interaction of a single atom with radiation field is a basic interaction in physics. In [6], a 
nonperturbative fully quantum-theoretical analysis describing the transient spontaneous 
emission of an initially excited two-level atom in a one-dimensional cavity with output 
coupling is presented. In [7], they present observations of quantum dynamics of an 
isolated neutral atom stored in a magneto-optical trap. 

The theoretic measure adapted in pE] is a typical one which takes the laser field 
as a classical field. However, considering the quantum nature of the driving field, one 
may get results which is different from the result under classical treatment. There are 
generally two ways to take account of the quantum nature of the field. One is to add 
quantum fluctuations to the classical treatment j8] . However, there are many operations 
in quantum computation, then in this case whether this method is suitable is still to be 
considered. The other way is to quantize the field and work out the result [9] . To do this, 
we first consider the Rabi oscillation driven by quantized pulse train. This is a basic 
atom-photon interaction process, and quantum computation is only one application 
of it. Then we can analyze the failure probability in ion-trap QC, and have further 
discussions. 

Rabi oscillation driven by a quantized continuous-wave (cw) field, accompanied by 
collapse- revival phenomenon [TQJ [Til 1121 EE], is a typical phenomenon of atom-photon 
systems. However, Rabi oscillation driven by quantized pulse train has not been fully 
investigated. Conceivably, it may have different phenomena compared with that driven 
by a cw field. 

Fault-tolerant quantum computation (FTQC) allows the computer to work 
normally even when its elementary components are imperfect. However, the threshold 
theorem in FTQC requires the failure probability of each components below some 
threshold [H]. Then we can compare the failure probability with the value of threshold, 
and reach some meaningful conclusion in quantum computation [9]. 

This paper is arranged as follows: in Section [2} we describe a method to deal with 
quantum transformation of a two-level system after one coherent pulse, which expresses 
the relationship between the density matrices for the two-level system before and after 
one coherent pulse. In Section [3] we investigate the properties of Rabi oscillation driven 
by pulse train. In Section H] we study this kind of Rabi oscillation in ion-trap quantum 
computation, and get the failure probability. In Section |5] we give out some discussions. 
In Section [6] some conclusions are reached. 
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2. Quantum transformation of a two-level system involving one coherent 
pulse 

2.1. Modeling 

The two-level system driven by repeated pulses is an open system, the usual way to deal 
with such a system is Kraus summation and master equation method. However, for the 
specific problem here, which cannot be easily solved with those methods, we study it in 
this way: after a single pulse, we get the density matrix for the whole system (including 
a two- level system and the laser field), then obtain the reduced density matrix for the 
two-level system. We can get the relation for the state of the two-level system before 
and after the pulse, then the state of the two-level system after repeated pulses can be 
get. 

In [U [15], they use the Jaynes-Cummings model (JCM) [16] for the case where an 
atom in the free space interacts with laser field. However, the JCM is the model to 
describe the interaction of atom and single-mode field in a cavity. Actually, there are 
some discussions [TTt [T8| IT9] on the validation of the JCM in the multi-mode case. For 
example, in the paper by Enk and Kimble [TH], at section 2.3 "Atom-light interaction", 
they considered the case where an atom in free space interacts with a laser field, making 
use of Hamiltonian of JCM in Eq.(lO). They also pointed out that the Hamiltonian in 
Eq. (10) in the paper is valid for atoms in free space for less than one Rabi period, but 
there is no strict proof. 

We analyze the issue as follows: the sources of decoherence can generally lead 
to certain amount of failure probability on a single qubit or a pair of qubits. After 
many operations on the same qubit (or the same pair of qubits), the failure probability 
will generally accumulate to reach the threshold in the threshold theorem of FTQC. 
The corresponding operation number is the upper bound of operation number in one 
error-correction period when the given source of decoherence exists. Each source of 
decoherence considered can generally give corresponding upper bound of operation 
number, but this upper bound can be enlarged by technique improvement. When the 
upper bound increased to a huge value, it is actually not a bound. 

The decoherence caused by field quantization can also give its upper bound of 
operation number. It is obtained theoretically, thus cannot be enlarged by technique 
improvement. The calculation of this decoherence should include the interaction of all 
modes in the radiation field with the two-level system. When using JCM, only one mode 
of the field is considered, and can also give an upper bound of operation number 0x- 
The accurate upper bound of operation number from field quantization /3 < flx, because 
the spontaneous emission induced by vacuum modes is not considered in JCM. Then if 
we use JCM to estimate the upper bound of operation number in one error-correction 
period from field quantization, we can get meaningful results. The two-level system 
driven by pulse train can be described as 

H = hg (e^a+a + e _i V<r_) , (1) 
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where g is the coupling constant, is the beam phase, o~ + and er_ are the raising and 
lowering operators of the two-level system, and a' and a the creation and annihilation 
operators of photons, respectively. Then the unitary time-evolution operation is given 
by 

U{t) = cos (gty/ala + l) |1>{1| + cos (gts/ata) |0)(0| 

sin [gtyfcfia + 1 ) sin ( gt\J a^a + 1 ) 

-a|l)(0| + e-'^t \ y |0)(l| 



V at a + 1 



(2) 



with |0) and |1) the ground and excited state of the two-level system respectively. 

Generally, the initial state of the whole system is \if) (0)) = J2^=o c n\ n ) ® 
(a|0) + where |c n | 2 = 6 " n , and \a\ 2 + \(3\ 2 = 1. A single qubit gate is usually 

implemented through a kir pulse in Cirac-Zoller scheme, whose duration t satisfies 
gt y/n~ = *y [T5], with n the mean number of photons in the pulse. After a kn pulse, 
the state for the two-level system and laser field is 




kir\/n.. , ;i .kitx/n^, 
( (*s(— l=r)\0,n) -ie^sin(— 1) 



2Jn 



+/?< J2 Cn 



. n=0 



COS 



\l,n) — le v sm( = — )|0,n+l) 



2y/h~ 



The corresponding density matrix for the state in (jHJ) is p^L; = 



• (3) 



This 



matrix contains the information of both the two-level system and the field, but we are 
interested only in the two-level system. Thus we obtain the reduced density matrix 

\a\ 2 S 4 + - a*(3)e^S 2 + |/3| 2 (1 - S 6 ) af3*S 5 + i(|a| 2 e i( ^i - \(3\ 2 e-'^S 7 ) + a*f3S 3 

a*pS 5 - i(H 2 e i( ^i - \/3\ 2 e-'^S 7 ) + a/3*S 3 \a\ 2 (l - S 4 ) - ±(a/3* - a*/3)e i *S 2 + \/3\ 2 S 6 



P 



(i) 



n=0 



n! V ti + 1 



n=0 

oo 

s 3 = E 



e n 



ft! 



n=0 



s 4 = E 

n=0 

oo 

S = = E 



e n 



n kix^fri krry/n + 1 , 
cos( — ) sm( 



2(n + r 



sin 



2^fn~ 2^fn~ 
,kn\/n + 1. 



n 



n n n / n .k^Jn .kny/n + 1. 

— I — \/ 7 sm ( 7^ J sm ( /= — 

n\ V n+ 1 2^ 2^ ' 



cos 



n! 



,kny/n 

~2~Vn 



n=0 



s. = E 



n=0 



n! 



k-nJn kny/n + 1, 

cos! =r- ) cos( = — 

v 2^n ' y 2^n ' 

. kiry/n + 1. 



cos 



2^n~ 



-)• 



(4) 
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n=0 



nl 



n .kiiy/n + 1 kiiJn 
~ cos ( — nr ) sm (^-7^)- 



2.2. Transforms of the density matrix after a coherent pulse 

Consider the relationship between p^ and the density matrix of corresponding initial 
state p(°) = \ip (0))(ip (0) \. For a two-level system, the density matrix p satisfies 
the condition p = ~(I + r • c) [2], r is the Bloch vector for state p, \r\< 1, 



a 



a x a y a z 



Let r( r 



(m) (m) (m) 



denotes the Bloch vector of p^ m \ An arbitrary 



trace-preserving quantum operation is equivalent to a map of the form r — > r' = 
Mr + c [H], here M and c contain the properties of the system and are independent 
of the state. Based on this, it can be seen that = Mr^ + c, here c = 



Sre-* - S 4 + S 6 -l 

M = 

then r( m ) = Mr {rn - 1] + c 



S 3 + S 5 
Mi 



Mi 



5 5 -5 3 -(e^i + e-^SV) 
5 2 e^ ^ 4 - 5 6 



Calculation of the sums in the density matrix 

It is necessary to get accurate values of Si (i = 1, • • • , 7) to evaluate the behavior of pulse 
train. The usual algorithm (saddle-point approximation) can only reach a precision of 
1 / Vn. Our algorithm achieving any given precision instead of the usual algorithm is as 
follows. 

Suppose n is not small, for the sum 
Y] , fio(n,n,k), 

n=0 

(1) Substitute n in fio(n, n, k) with (x+l)n , we get fn(x, n, k) = f i0 ^(x+l)n, n, kj . 

(2) Do Taylor expansion to x p for fn(x, n, k) at x = 0, and get /^(x, n, fc). 

(3) Since sum 5^^1 6 """ n fc can be get accurately, we replace x in fi2(x,n, k) by 
2=2 and get f a (n,n,k). 

(4) Use f i3 (n,n,k) instead of fio(n,n,k) in the expression of Si(n,k) to calculate 
the new sum and get /^(n, k). 

(5) Substituting n into fu(n,k), we obtain a high-precision result of the original 
sum Si(n, k). The value for Si (i — 1, • • • , 7) in the cases where we expand fa(x, n, k) 
to a; 10 and x 15 are compared in Table HJ 

The precision of the sums Si(i = 1, 2, • • • , 7) is ensured by the following theorem: 
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Table 1: Values for Si(i = 1, 2, • • • , 7) for n = 10 4 and k = 2. "Valuer" denotes value 
of the resulting sums of the algorithm when we expand fn (x, n, k) to x 10 and "Value2" 
denotes that when we expand fa(x, n, k) to x 15 . Valuel and Value2 are the same to the 
precision 10 -23 . 



Sum 


Valuel 


Valuc2 


Si 


0.000 039 303 916 656 063 668 561 519 091 


0.000 039 303 916 656 063 668 561 194 770 


s 2 


0.000 039 265 164 255 300 772 996 074 590 


0.000 039 265 164 255 300 772 995 750 283 


s 3 


0.000 246 659 192 761 352 167 541 307 293 


0.000 246 659 192 761 352 167 542 402 758 


s 4 


0.999 753 309 972 685 637 856 777 333 369 


0.999 753 309 972 685 637 856 776 237 858 


s 5 


0.999 753 316 133 881 571 308 212 070 145 


0.999 753 316 133 881 571 308 210 974 684 


s 6 


0.999 753 322 301 165 250 291 025 614 276 


0.999 753 322 301 165 250 291 024 518 866 


s 7 


0.000 039 226 416 698 193 975 826 600 887 


0.000 039 226 416 698 193 975 830 095 264 



Theorem 1: For every given integer I « n, let 
In (y2n l -^(l + l)lnn) 



P 



«o 



\ In n — In ( (7 + 1) In n 



(5) 



n 



(, + i )M + j(i ± fH !+2(1 + 1)hi 



n 



If a < a « \/n, then 

yv e — _^ _ e — n _^ + / 1 

n=0 n=0 

here p, fio(n,n, k), fi^(n,n, k) are parameters defined in the algorithm above. 



(6) 



Theorem 1 can be proved using the following three lemmas (see Appendix A for 
the detailed proof): 

Lemma 1: For every given a « sfft, 



E 



711 



-{fio{n,n, k) - /i 3 (n,n, k)) = o 



a 



P+i 



(Vn) 



p-i 



(7) 



n=n— ay/n 

here p, f i0 (n,n, k), f i3 (n,n, k) are parameters defined in the algorithm above. 



Lemma 2: For every given integer I « n, a « \/n, if a > v/^T+T) Inn, then 

1 



E 

n=0 



e -r < —i 

nl rv 



(8) 



where k = [~n + o;\/n] . 

Lemma 3: For every given integer I « n, a « y/E, if 



1 (Z + l)lnn /(/ + l) 2 (lnn) 2 nn iN1 

n > ^= + - j= + \ — ^ + 2(/ + l)lnn, 

'n Vu V n 



(9) 



then 
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where k' — [n + oja/^J . 

We expand f\(x,n, k) to x 15 at x — 0, and find the value of the sum is the same 
to the precision 10~ 23 as that when we expand fi(x,n,k) to x 10 . It implies that the 
precision of the sum is much higher than Eq. ([5]) shows. The reason is probably that we 
have not considered of the periodicity of trigonometric functions, and the precision of 
the sum may be remarkably improved since the positive and negative terms will cancel 
each other out. 

For small n, we need only to require t satisfying 

it- 1)! > e- fl n t+l , (11) 

where t is the parameter in sum Si(n, k) = Yln=o £_ ^f~^°( n ' ^> ^) + °\Wj ■ ^ or a & Yen 
precision I, we can search out the smallest t satisfying (11 II) . e.g., when n = 10 and 
I = 20, we set t = 55. 



3. Population inversion 



3.1. Final state of the two-level system after pulse train 
Provided = Mr {m ~ 1] + c , then 

r (m) = M m r (0) + ( M m-l _| + M + I)c. 

It can be seen from Sec. 12.21 that 



M r - 



For any real matrix A 



(S 3 + S 5 ) m 
O M™ 



S 5 -S 3 -(St + Sr) 

S2 5*4 + Sq — 1 



a b 
c d 



we obtain (see Appendix B) 



Am ^ r A - 



(m) 



2ig 



-X 26 
2c if 



where A 



(m) 



(12) 



(13) 



± 



ym _|_ ym w }th \ l anc [ ^ 2 eigenvalues of A, K = d — a,Q 



-i\/(a - d) 2 + 46c. When (a - d) 2 + Abe < (which is the case for Mi) 

J 



A m = \X\ m cos(m6)I + sin(m6 



VdetJ 



where |A| 2 = ad - be, sin9 = \ J2 - a2 +£+f c , J 



m—l 



a — d 2b 
2c d — a 



. Therefore, 



m—l 



I + Mi + ■ ■ ■ + M™- 1 = \J2\X\ j cos(j9) I + [ \ X \ j sm(j^) 



J 



j=0 



j=0 



VdetJ 
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m— 1 



m—1 



m—1 



we have 



3=0 3=0 3=0 

1 - |A|e- ie - \\\ m e< m °) + \\\m+i e i( m -i)e 
~ 1 + |A| 2 — 2|A| COS6* ' 

^ . _ 1 - |A| cos^ - \X\ m cos(me) + |A| m+1 cos(m- 1)6 

^ |A| COS( ^ } " l + |A| 2 -2|A|cos# ' 

. , _ |A| sin<9 — \\\ m sm(mO) + |A| m+1 sin(m- 1)6 
? |A| Sm( ^ } = l + |A| 2 -2|A|cosfl ' 



thus 



I + Mi + • • • + M[ 
1 



m—1 



(1 - |A|cosfl 



1 + |A| 2 — 2|A| cos6» 
+ (|A|sin0- |A| m sin(m0) + |A| m+1 sin(m - 1)6) 

= B[ m) I + Bi m) J, 



cos(m6) + |A| m+1 cos(m - 1)6)1 
J 



VdetJ 



then 



r O) — 



+ 



(s 3 + s 5 y 
o 

l-(S 3 +S 5 ) m 
1-S 3 -S B 



cos(m6)I + sm(m6)^=L= 



,(o) 



O B^I + B^J 



c. 



3.2. Population inversion after pulse train 

Suppose the initial state is |1), if we have applied kn pulses for m times, the population 
inversion is 

Wm = 1 (1 _ r M) _ 1 (1 + r M) = _ r (m) > 



we have 
W m = |A 



cos(m6) + sin(m#)- 



322 



v^etjJ 1+ |A| 2 -2|A|cos£ 



x{ [|A| sind - |A| m sin(m0) + |Ar +1 sin(m - l)d\ / 21 (S 7 - S ± ) 
L J V det J 

+ [(l - |A|cos0- |A| m cos(m0) + |A| m+1 cos(m - 1)^ 



+ (|A| sind - |A| m sin(m0) + |A| m+1 sin(m - l)fl) -jt= = ] (S 4 - 5 6 )}. (14) 
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To get the inversion between the mth and (m + l)th /c7r pulse, we should first get 
the corresponding density matrix of the two-level system 

Pm {t) = ti{U{t)[p^ ® Pl ]U\t)}, (15) 

where U(t) is the unitary time-evolution operator mentioned earlier, and p\ is the density 
matrix for the laser field. A detailed calculation gives the probability that the ion is in 
state |0): 

where S, ^ = £~=o e -^f cos 2 gty/E, S 9 = £~ sin 2 gty/K+T and S l0 = 
Yl^=o 6 rX y/n s ^ n ^Qfcy/n, which can all be got with high precision using our algorithm 
in Section [231 and the inversion is W m (t) = 1 — 2p(t), < t < 

Now think of the difference of the oscillation driven by pulse train and a cw field. 
The population inversion for repeated 2tt pulses is shown in Figure [la] (given n = 10). 
We find there is a dual-pulse structure in every period, where the amplitude starts to 
increase from the point of 2tt pulses. The inversion decreases exponentially, different 
from a Gaussian function collapse envelope driven by a cw field. Besides, there is no 
revival phenomenon, but a small nonzero amplitude exists. The reason for this behavior 
is probably that the system we consider here is an open quantum system with dissipation. 

The inversion at the points of 2ir pulses when n = 10 4 is plotted in Figure El 
Results of fitting is 1.0031e- a0002ArR , 1.0193e- a0003ArR , 1.025e- a0005ArR for k = 1/2,1,2 
respectively, here is the number of Rabi periods. 



Pit) = « (Ss + S 9 ) + ri m \t)(S 8 - S 9 ) + r™(t)S w 



4. Failure probability of gate operation realized through Rabi oscillation 
driven by repeated pulses 

4-1. Estimation of n 

n which determines Si is an important parameter in our discussion. To estimate the 

mean number of photons in one pulse, we assume a fictious pulse propagating in 

the opposite direction at the same time. They form a standing wave instantly while 

overlapping in space. It can be seen that the mean number of photons in each pulse is 

about half of that in the standing wave. We now focus on mean number of photons in 

the imagined standing wave. 

The electric field E can be expressed as E = £y/h~. £ is usually given as £=\J ^7 

[20] . where u is the frequency of the single mode in a cavity, and V is the volume of the 

cavity. It can be seen that V ~ Act, with A the cross-sectional area of the beam, thus 

n = ^j^E 2 . For a kir pulse, gt\ffi = g ~ j- = -0=, with p ~ eao the electric dipole 

moment of the ion, e the charge of a electron, and Bohr radius, then we get t = 

Then we have n = ^^E. 
, p 

When the laser beam interacts with an ion, not all the photon in the beam are 
effective. Think of the effective mean number of photons. When a laser beam is 
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Figure 1: Population inversion driven by different fields, given n — 10, r = gt. (a) 27r 
pulse train case. There is a dual-pulse structure in every period, where the amplitude 
starts to increase from the point of 2ir pulses. The inversion decreases exponentially, 
different from a Gaussian function collapse envelope driven by a cw field. Besides, there 
is no revival phenomenon, but a small nonzero amplitude exists, (b) Corresponding 
inversion driven by a cw field. 

applied to a trapped ion, the total resonant scattering cross section for an atomic dipole 
transition is a = 3A 2 /27r |21j . and the cross section for scattering out of the paraxial 
modes is a e ff = 3\ 2 /8n [22]. Then the effective interaction area is & e ff, and the photons 
in volume cr e ffct is effective. For each photon, the probability to be in area cr e // is 
and the probability are independent for the photons. It can be seen that the effective 
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k=1/2,l=4 

— k=1,l=4 
k=2,l=4 




0.2 



1000 2000 3000 4000 5000 6000 7000 

Figure 2: Inversion at the points of 2tc pulses E(W) versus number of Rabi periods 
N R , n = 10 4 . Fitting results are 1.0031e- a0002AfR , 1.0193e-°- 0003iVR , 1.025e- a0005JVR for 
k = 1/2, 1, 2 respectively. 



mean number of photons is 

_Cr e ff ktO&eff^rn 

n e ff = n — J -± = 1J — E. (16) 

A A p ( J 

One case people are interested in is the sideband transition, where the laser detuning 

A = ±u; t , here u t is the frequency of the trap. Because of AC-Stark shift and off- resonant 

transitions, the sideband Rabi frequency Q + has upper bound [23]. People have adopted 

methods to partially cancel the effect, and it seems feasible to have Q + < u t for special 



temporal and spectral arrangements of the laser field [21]. Since fi + = jy 2MujS *"> w ^ n 
M the mass for a single ion, we have 



A 2M 3 



From [25] and [26], it can be seen that 
ea$E pE 



Q - 



e 2 



Ane Mz^ ' 

where z s is the order of the separation between ions and is typically 10 to 100 /iin. 
Suppose z s = £A, from ( TITl) and ( fT8l) . we get 



2y/2h, e 



2 



5 



E<— — ( )?M^f?A"5. (19) 

pn Aneo 



Substitute back to ( fTBI) . we get 



32a2vrx V e 

= 6 x 10 7 fcM-^-fA3. (20) 
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In the cases we consider, it is suitable to limit k < 2, 9u < M < 200u (u= 
1.66057 x 1(T 27 kg). For M = 9u, k = 2, we get 

n = 3.4 x 10 14 r^i 

We can see that a large A and a small £ result in a large n. The curves of lg(n) is plotted 
in Figure |3] versus parameter £ from 2 to 100. When A = 10~ 6 m and £ = 2, we get 
n = 2.3 x 10 3 . 




Figure 3: Logarithm of mean number of photons lg(n) as a function of £ and A. It can 
be seen that n increases with A and decreases with £. 

There are also authors calculating n in a kir pulse in another way [15] : they 
considered the situation where a single laser is used to drive Rabi oscillation of the 
atom, and adopted the formalism introduced by Blow et al. [27], taking the laser 
as a continuous-mode coherent state. They worked out the interaction time t for kir 
pulse as t = ^r^, where d is the coupling constant of the atom and laser, and 

P is the power of the laser. Thus, the mean number of photons in one kn pulse is 

™ ~ /£z = utji\J wnere is the frequency of the representative single-mode 
coherent state. Thus, obviously, they take all the photons in area A as effective photons 
when considering the interaction, but actually the number of effective photons is much 
smaller. 

4-2. Accuracy of gate operation 

Suppose we have applied m times of coherent pulses and reached a state = 
i(7 + r^" 1 ) ■ a). Let |^) = a|0) + /3|1) be the expected state, the accuracy rate of 
gate operation realized through Rabi oscillation is 

p^ = ( a *(0\+P*(l\)p( m \a\0) + (3\l)) 

= i (\ _l r (o) r M + r (o) r M , (0) {m)\ 
o \ y y J 

= i(l + r (o). r M) ; (21) 
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for a mixed state, |r^| < 1, then p s < 1. The failure probability is 



Pi 



P 



13 

M A 



detailed calculation results (see Appendix C ) 



Pi 



(m) 



-(r 



(0) . r (m) 



(r^) 2 [(S 3 + S 5 ) m - 1] + (H 0) ? + (rl 0) ) 2 )[|A|^ cosM) - 1] 

+ |A|^ sin(m#)(det J)"* [((rf) 2 - (rf) 2 )(a - d) + rfVf (b + c)\ 

+ i?i m) (rj 0) c y + rfc z ) + Bi m) [(r^c y - rfc z )(a - d) + rfV + ^ 0) c z 6] 

Then we average over all initial states of the ion, and get the average failure 
probability. The failure probability for kir pulses with different n are shown in Figure HJ 
It can be seen that the failure probability increases with the number of Rabi periods 
N R and the value k, and is inversely proportional to n. 



(a) 



(b) 



k=2,l=4 

- k=1,l=4 
k=1/2,l=4 




100 




80 100 



Figure 4: Failure probability pt versus number of Rabi periods A^r, kn pulses are applied, 
and n is 10 4 (10 6 ) in a (b). The failure probability pf increases with the number of Rabi 
periods Nr and the value k, and is inversely proportional to n. 



5. Discussions 

5.1. The reason for tracing out the laser field after each pulse 

For the Cirac-Zoller CNOT scheme, there are five steps to implement it. The laser 
phase in the five steps are 7r/2, 0, 0, and — vr/2, which are correlated. However, there 
is no correlation between two different CNOT gates. Then after completing one CNOT 
gate, it is valid to trace out the laser field. The calculation in Sec. |2]for different pulses 
is an analog of different CNOT gates. This treatment is meaningful for our purpose, 
i.e., to estimate the lower bound of failure probability, or the upper bound of operation 
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number, after many CNOT operations. More-detailed analysis of the lower bound of 
failure probability after Cirac-Zoller CNOT operations can be seen in [28], which is to 
be published. 

5.2. The permitted depth of quantum logical operation 

Failure probability we have calculated for the ir sideband transition amounts to 10 -2 
after approximately 10 2 operations when n = 10 4 . For controlled-NOT (CNOT) gates, 
there are five steps in Cirac-Zoller scheme, and two steps are realized via Rabi oscillations 
driven by n pulses. Generally speaking, the failure probability after five steps is not 
less than that after one ir pulse. Then the failure probability after repeated 7r pulses 
is a lower bound of the failure probability after repeated Cirac-Zoller's CNOT gate. 
Thus the lower bound of the failure probability is 10 -2 after approximately 10 2 CNOT 
operations when n = 10 4 . 

The threshold theorem in quantum computation declares that an arbitrarily long 
quantum computation can be performed reliably if the failure probability of each 
quantum gate is less than a critical value. Knill used numerical calculations and obtained 
a threshold of the order 10~ 2 [29] based on a fault-tolerant structure suggested by 
himself. P. Aliferis et al. reached a threshold of 1CT 3 with provable constructions [30] . 

A parameter called permitted depth of logical operation describing the property 
of a physical realization scheme of QC has been given [9]: considering that different 
number state components of the driving field lead to different oscillation amplitudes, 
which become uncorrelated gradually, we can see that failure probability of quantum 
logic gates has a theoretical limitation. Combing this limitation given by the quantum 
nature of the field with the threshold theorem in FTQC, we can get permitted depth 
of logical operation. This parameter limits the number of operations on any physical 
qubit in one error-correction period. Then the permitted depth of logical operation here 
is less than 10 2 . 

5. 3. Others ' proposals which may have different results 

For Rabi oscillation driven by microwaves, the failure probability may be much smaller 
because of a large mean number of photons, but it becomes difficult to individually 
address each ions. Although an additional magnetic field gradient applied to an 
electrodynamic trap may individually shift ionic qubit resonances [31], thus making 
them distinguishable in frequency space, whether it can improve the permitted depth 
of logical operation needs further investigation. 

There exists a two-qubit gate scheme totally different from the Cirac-Zoller gate, 
such as the scheme implemented by the NIST group [32]. In the scheme off- resonant 
excitations of the stronger carrier transition is absent, and this allows a greater gate 
speed thus a higher laser intensity. Besides, additional Stark shifts can be efficiently 
suppressed by choosing almost perpendicular and linear polarizations for the laser beams 
[33] . Hence, studies on this type of gate may lead to different results. 
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6. Conclusions 



Firstly, we investigate Rabi oscillation of a two-level system driven by pulse train. We 
develop an algorithm to solve the infinite summation, with a high precision which has 
never been reached by former authors. We find that in this kind of Rabi oscillation, there 
is a dual-pulse structure in every period. The envelope of population inversion collapses 
exponentially, different from a Gaussian function collapse envelope driven by a cw field. 
Besides, there is no revival phenomenon, but an approaching to a tiny platform. 

Secondly, we consider the application to gate operation in ion trap quantum 
computation. We give a lower bound of failure probability. Our result is: that when 
the wavelength of the driving field is of the order 10~ 6 m, the mean number of photons 
cannot be greater than 10 4 . Then, after about 10 2 CNOT gates in Cirac-Zoller scheme, 
the lower bound of failure probability is of the order 10 -2 . 
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Appendix A. PROOF OF PRECISION OF THE ALGORITHM IN SEC. 

Proof of Lemma 1: For n — a\M < n < n + a>\M, i.e. — % < x < -7=, after 
expanding fn(x, n, k) at x — 0, we get the result f i2 (x,n,k) satisfying f i2 (x,n,k) = 
fio(n,n,k) + o(x p ). It can be seen that fi 2 (x,n,k) = f i3 (n,n,k), thus we have 
fio(n, n, k) = fz{n, n, k) + o(x p ), then Eq. (j7j) is proved □ 

Proof of Lemma 2: For every given n satisfying n < k + 1 < n, we have 

< (fc+l)!> tilUS 



k _„• k 



f.k+1 ^k+l 



n=0 n=0 v ' 

k\k. 



From Stirling's formula k\ = y2itk(~) ei2fc,0 < 6 < 1, we have 

e ~"^TT < (l) k e- n n k+1 = e fc ""n(-)- fc . (A.2) 
Substitute k in formula ( 1A.2I) with n — a\/E, we have 



n Jn 



When a > y^T+T) Inn, we have ne ° 2 < ^, inequality ()8]) is proved. □ 
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Proof of Lemma 3: It can be seen that 

00 ^k' 00 ^n-k'lJ\ 

e ~^\ =e Tn S 



n\ k'\ t-^, n\ 

n=k' n=k' 



n k , = n n k ' kf + 1 
k'\ ^ { k' + V k'\ k' + l-n 

n=k' 



when k' > n + 4, i.e, , , fc t 1 _ < k' we have 



~-j /C ~)~ 1 . ; .V I A/ 1 



e "— — < — e 'V Ti <ne 



fc'! v A;' y (fc'-l)! 
with Stirling's formula we get 

( e )fc ' e -n-fc'+i < - e -n ™* < ne ~n n k'-i = ne -n n k'- K 

Let A = < 1, 77 = n^r, we then have 

ne- n n k '- 1 = ne^n*'- 1 < (t^)^ 1 < 4tt ^ (")* - eA > 

rr k — 1 n l+l i] 

^ A(l - In 77) > 1 + In A. 

Let A = 1-A, withO < A < 1, from ln(l+x) < rr-irr 2 (x < 0), we get In A < -A-|A 2 , 
then a sufficient condition of (^) A — eA > is: 

(l-A)(l-lnr ? )>l-A-iA 2 , 

which results in A > A , here A = — In 77 + \j (In??) 2 + 2 ln?y. Let = n + ckoV^, 
we get 

1 



«o = — 7E 



1-A 







- 1 )n+ 1 



1 / 1 + 1 , /// + 1 , _\ 2 J + l 



— + I — = — In n + \ / ( — = — In n ) + 2 — = — In n | \/n 
Vn \ ™ \ \ n J n 

by using = 1 + In rj + -y/(ln w) 2 + 2 In rj and i] = (n)~*~ . Because A > A a > a , 
we get a sufficient condition of Lemma 3: 




1 (Z + l)lnra /(/ + l) 2 (lnn) 2 n/1 

— + - £ + \ — - + 2(1 + 1) Inn. 

n vn V n 



then Lemma 3 follows. □ 

Proof of Theorem 1: From Lemma 1,2 and 3 we get: for every given I « n, 
if a satisfies 



Full quantum treatment of Rabi oscillation driven by pulse train 



17 
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Appendix B. CALCULATION OF Mi m 

Let 



where 



then 
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Figure Bl: A(r) = (a — d) 2 +4bc < versus r = gt, where t is the pulse width. Different 
A(t) results in different behavior of Rabi oscillation driven by pulse train. For the cases 
we consider, r < 1, we can see A(r) < 0. 



Let (1, X2i) T and (1, X22) T be the eigenvectors of Mi with corresponding eigenvalues 
Ai and A2, then 



X21 = —[d — a + \JJa — d) 2 + 46c] , 



2b 



X22 



hd-a- v / (a-rf) 2 + 46c], 



2b 



Ai = -[a + rf+ \/(a - d) 2 + 46c], 
A2 = ~[a + - A/( a - + 46c], 



we have plot A(r) = (a — d) 2 + 46c versus r = in Figure EH A(r) is below zero for 
the cases we are interested in (r < 1). 
Let 
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Denote A™ ± A™ = Ai m) , d-a = K, ^{a - d) 2 + Abe = iQ, we have 
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It can be seen that |AjJ = |A 2 |, let Ai = |A|e ie , A 2 = |A|e _ie , using |A| 2 = AiA 2 = ad — be 
we get 

A^ m) = 2(ad-bc) m/2 cos(m#), 
A (m) = 2i(ad - bc) m/2 sin(m0), 

where 9 satisfies sin 9 = 



2ad-Abc-a 2 -d 2 
4(ad-bc) ' 



then 



M™ = |A| m cos(m9)I + sin(m6 
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VdetJ 



where 



J 



a — d 2b 
2c d — a 



(B.2) 
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It can be seen from (fI2]h (fill), flRTD . flR2|) that 
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